Of Bombs and Boats and Mice and Men 
A random tour through some scaling laws 

Niall MacKay 

I. I'll begin with a true story which has become semi-mythical. Shortly after World War Two, the 
US Atomic Energy Commission released a film of the 1945 Trinity atomic bomb test. The energy 
yield remained secret, having been estimated only with some difficulty. So the Americans were 
most surprised when the British fluid dynamicist Geoffrey Taylor published, in 1950, an accurate 
estimate merely by studying the AEC's pictures [1]. Embellishments of the tale have the Americans 
wondering how on Earth he did it, with the CIA visiting his house in the middle of the night to 
search through his papers. 




Figure 1: the Trinity atomic bomb test. 

The way the story is often told nowadays, his technique was simple dimensional analysis. Suppose 
that, at time t, the radius r of the blast wave depends only on t (the time elapsed since the 
explosion), on the air density p, and on the energy E released in the blast. (That is, once the blast 
wave is propagating it knows nothing of the nature of the explosion which caused it.) Energy is a 
mass multiplied by the square of a speed, which we write [E] = ML'^T^'^, while for the density 
[p] = ML~^. So to cancel the masses, we must combine E and p as E/p. This has dimensions 
[E/p] = L^T^"^, so that to form a length from this and t we must take (Et'^/p)^/^, and finally 



/Et 



where C is a dimensionless constant, universal for such waves, which can be estimated from 
conventional explosive blasts. Used with pictures such as that in Fig. 1, which helpfully gives 
length and time scales, this formula enabled Taylor to estimate E. 

This calculation is striking especially for the unexpected appearance of the one-fifth power — after 
all, when we clap our hands we are used to the effects propagating very differently, with a fixed 'speed 
of sound' V, and r oc vt. It is often used to illustrate the power of dimensional analysis, but this 
is really another myth, potentially more pernicious than the more exciting embellishments, for the 
dimensional result only emerged as part of a deeper analysis which included both the mathematics 



of self-similar solutions of the relevant partial differential equations (PDEs) and a careful treatment 
of the physics involved. The dimensional derivation is true, but hardly independent^ 



Another classic example is the problem of roasting times. How much longer does it take to roast 
a 5kg turkey than a 1kg pheasant? Most recipes give a linear formula: so many minutes per unit of 
mass, plus a fixed time. However, this is only an approximation to a scaling which can be discovered 
by dimensional analysis. Let us assume that the birds are similar in shape, roasted at the same 
temperature, with similar flesh, and that we wish their centres to reach the same temperature to 
be properly cooked. Then the cooking time t depends only on the typical distance I from surface 
to centre and on one physical parameter, the thermal diffusivity k. This is associated with the one 
microscopic physical law which we need to solve the problem, the 'heat equation': that the rate of 
change of temperature is proportional (via k) to its Laplacian (of second spatial derivatives), so that 
[k] = L^T^^. (We don't need to consider precisely how n depends on the thermal conductivity, 
density and specific heat capacity of the flesh.) This immediately gives us that t ~ l"^ /k. Since the 
bird's mass m ~ l^, we have t = Crr?!"^ for some dimensionless constant C . So, if the roasting 
time for the pheasant is ~ Ihr, then 



Dimensional analysis, then, is fundamentally bound up with the concept of scaling, which is 
crucial in applying mathematics. For example, suppose we wish to scale up a laboratory chemical 
experiment to an industrial process. If the reaction produces heat then we see immediately that 
the apparatus may not simply be self-similarly scaled up, for then the extra heat produced (in 
proportion to the volume of reagents) will not be matched by the ability to dissipate it safely (in 
proportion to the surface area of the vessel). Rather we should identify the dimensionless numbers 
which characterize the problem and keep these invariant under the scaling. In the examples above, 
there is only one dimensionless number, Et'^/pr^ for the bomb blast and nt/l'^ for the roasting 
times. If there were more, we could write down in scale-free form the equations which relate them. 

To keep a dimensionless number constant is precisely to impose a power-law scaling, as in ([T]). 
Indeed for scale-invariance r must scale as a simple power of each of E, p and t, for this is the 
only functional form which is invariant under changes of scale (or, equivalently, of units). With 
r = Kt^/^ (taking K = C{E/ p^/^ to be constant), multiplying t and r by constant factors merely 
changes K; it does not change the form of the equation, the power law. It was precisely this 
similarity property of the blast wave that enabled Taylor to simplify the solution of his PDE. 

The basic picture we have described holds in all applied mathematics, and dimensional analysis is 

''in fact Taylor had been working on the mathematics of blast waves throughout the war, publishing his results in a 
special, classified volume of the Proceedings of the Royal Society. Thanks to Martin Smith for drawing my attention 
to this volume. 

^We read y ^ x'^ as 'y scales as x^' , and mean (with no more precision than is necessary for our purposes) that 
yx~^ is approximately constant. If the reader baulks at the use of I and the assertion that m ~ , we note that one 
could instead use rn and density p, with [p\ = ML~^ . 




and the turkey's roasting time is 5^/'^ ~ 3 hours 
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an essential tool, not only for all natural scientists but for anyone who applies nnathematics. Even 
money, for example, can be considered a primitive dimensionful quantity, with its own units (of 
currency). Yet the subject is taught only patchily to users of maths: in English schools it appears 
(as far as I am aware) only in the more advanced A-level mechanics modules. There is a particular 
danger, as we shall see in Sect. IV, in the practice of not including the units in the algebraic 
variables — that is, for example, in writing a speed as v ms^^ with v a (dimensionless) number, 
rather than simply as v, a dimensionful quantity with no intrinsic number but which acquires one 
when divided by some suitable unit vq, such 3S vq = 1 ms^^. It is striking that this is practised in 
the same examinations which test elementary dimensional analysis — surely a recipe for confusion. 

II. 




Figure 2: a 1' twig and a 100' tree. 



An impression is sometimes given in popular science that scale-invariance is most naturally realized 
as a property of branching networks — so that a twig, for example, has much the same architecture 
as a tree (Fig. 2). But of course we can find it in any process which is invariant under rescaling of 
a dimensionful quantity. Restricting for the moment to lengths, how about the pictures in Fig. 3? 
With modern materials and handling gear, even a large yacht can have a single large mainsail, and 
such 'Bermudan sloop' design is approximately invariant over an order of magnitude. 




Figure 3: a 14' dinghy and a 109' yacht. 



So what power laws might we seek for yachts? An obvious one is speed: how does the maximum 
speed of a yacht scale with its length? Well, /7the speed depends only on the boat's length I and on 
the acceleration g due to gravity, then the only way to extract a speed from these is as the square 



3 



root of their product: v ~ y/gl. This, of course, is where the mystery of the game of dimensional 
analysis lies: which variables should be invited to the party? Why g, but no property of the hull 
materials or of water? In this case the answer is straightforward: every 'displacement' hull, which 
makes progress by moving through the water, sets up a bow wave, of wavelength proportional to 
the boat's waterline length, and which it cannot overtake. A surface wave of wavelength A on deep 
water has speed y^gX/lir, and the result follows^] But such questions are generally very subtle [3j. 

So let's be flippant instead. How fast does a yacht's price scale with length? There's no underlying 
physical law here, just the kind of question to which a mathematician at a boat show, astonished 
by the high prices, wants an answer. 




log( 1 / ft ) log( 1 / ft ) 

Figure 4: how the price p of secondhand yachts scales with their length I: 
(a) raw data p ~ /3.8±o.2^ ^R^ = 0.71; (b) age-adjusted p ~ /3.5±o.ig-o.03(age/yr)^ = 0.87. 

The simplest way to look for a power law y = Cx^ is to take its logarithm^ so that logy = 
/31ogx + logC. Rescalings of x and y are translations of the origin. A power law is seen as 
a straight line, and the scale-invariance property is that a straight line remains straight under 
changes of origin. Fig. 4(a) shows the prices of some used sailing yachts [4], and (because they 
were secondhand) Fig. 4(b) removes the best-fitting exponential decay of price with age, about 3% 
of value per year. Both the high resulting power (price p ~ l^-^) and the quality of the fit (nearly 
90% of variance) are striking. But notice that, for a given length and age, there's still typically an 
e-fold variation in prices — due not only to condition, design and desirability, but probably also to 
the enormous cost of the kit (rig, engine, electronics etc.) on board. 

III. To an applied mathematician, the Body Mass Index (BMI), your mass m divided by the square 
of your height h, is a peculiar construct. It's clearly not invariant under self-similar scaling — why 
not the cube of height? — but then perhaps neither is a person's health: a six-foot man with the 

^For a classic introduction see [2]. The wave speed does not depend on properties of the water, and is typically 
valid for any incompressible, inviscid fluid. (Note that viscosity would introduce another classic dimensionless quantity, 
the Reynolds number.) The scaling of ship speed is quite accurate: a 25' yacht has a maximum speed of about 6 
knots, and a 600' aircraft carrier, however powerful, has a maximum speed of about 30 knots. 

''We use natural logarithms, but any base will do. 
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body shape of a toddler is probably not in good shape. In fact the reason for the square is empirical: 
among adults, levels of fat are maximally correlated with mh^^ when /3 ~ 2 [5]. But for children, 
where there's some genuine scaling going on, it seems highly inappropriate to use BMI, whose 
meaning will then be size-dependent. So, if m ~ /i^, what is (3 for children? 

An old argument for quadrupeds regards the animal as a beam supported by its legs, and derives 
the scaling law which places an upper bound on the proportionate sag of the beam [6]. If we 
apply this instead to a human, regarding the pelvis as a supporting beam for the torso, one finds 
7/3 < P < 8/3, probably at the upper end of the range [7]. For boys aged 5-18 this matches the 
age-binned data rather well (Figure 5)§| 



log( h / metre ) 
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Figure 5: how mass m scales with height h for boys in 
(a) UK: m ~ /i2-70±o.05^ Hong Kong: m ~ /i2-66±o.05. ^oth ^ > 0-995. 

IV. It's remarkable how well the basic structure of mammals scales: the dormouse and the bear in 
Fig. 6 share the same basic architecture, even though they are four orders of magnitude apart in 
mass. So what scaling laws might apply to mammals generally? 





Figure 6: a dormouse (mass ~ 20g) and a bear (200kg). 



^Growing girls have a power closer to three, due to their optimization for reproduction. 
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In his classic essay 'On being the right size', J. B. S. Haldane noted that 'You can drop a mouse 
down a thousand-yard mine shaft; and, on arriving at the bottom, it gets a slight shock and walks 
away, provided that the ground is fairly soft' [8]. For a first treatment of this, let's presuppose 
self-similar scalingl§ so that for an animal of length scale / we have m ~ . Air resistance is 
proportional to cross-sectional surface area (and thus to f') and to the square of speed v. Then, at 
terminal velocity, with air resistance balancing gravity, we have t^v"^ ~ and thus v ~ ~ m 
If the bear's terminal velocity is roughly 150mph, then that of the mouse is roughly 30mph. To do 
the analysis properly we should further consider the scaling of the process of sudden deceleration, 
but the effect is the same: as Haldane put it, 'A rat is killed, a man is broken, a horse splashes'. 

A classic scaling law for animals is that of metabolic rate s (the animal's power consumption) with 
mass m. A simple self-similarity argument, with the animal losing heat in proportion to its surface 
area, gives s ~ m^/^, but empirical data suggest s ~ m^/'^, 'Kleiber's Law' [lOj. With modern 
data, of carefully measured basal metabolic rates, and for marsupials (which are less variable in 
their exploitation of new energy sources than other mammals), the fit is amazing (Fig. 7(a)) [llj. 
One persuasive explanation [12] derives Kleiber's law from the scaling of branching networks, which 
nature uses widely: think of alveoli, vascular networks, tracheoles, trees. Lungs, for example, are 
effectively neither three- nor wholly two-dimensional, rather filling the lung cavities with as much 
oxygen-grabbing potential as a branching network can manage. But (returning for a moment to 
flippancy) we could instead create a connection with Sect. Ill by noting that if mammals' scaling is 
not self-similar but rather m ~ l^^'^, but their surface area is still dominated by terms proportional 
to P, then the heat-loss-through-surface argument immediately gives s ~ ~ m^/^ — Kleiber's 
Law! 
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Figure 7: how metabolic rate s scales with mass m for 
(a) Marsupials: s ~ mO-^^io.oi^ ^i?2^0.99; (b) Eutheria: s ^ ,7iO-72±o.oi^ 



0.96. 



This is a good point at which to answer a fundamental question about data fitting: are we correct 



^Clearly the eyes [9J and legs [8J scale differently from the torso. Biologists refer to non-self-similar scaling as 
'allometry', in contrast to 'isometry'. 
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to fit a line to the logarithmic data, rather than a power to the untransformed, natural data? The 
question is really about how we weight the residuals — in the extreme cases, about how we treat 
the outliers. Suppose there were a species of mouse with ten times the metabolic rate expected for 
its size, and a species of bear with a rate 10% greater than expected. Which is the outlier? If we 
believe that there is a scaling relation, a similarity in the architecture of creatures of such disparate 
size — and clearly we must do so if we're fitting a power — then it's the mouse that's the outlier, 
the stranger animal. This is treated correctly by the logarithmic data, where the mouse is about 
24 times further from the fitted line than the bear. But in the natural data it's the bear that's the 
outlier, about 16 times further from the curve|l| 

For the eutheria — the 'good beasts', the non-marsupial mammals — the picture is not quite so 
clear (Fig. 7(b)), and indeed the subject of metabolic scaling has become surprisingly controversial 
in recent years |[13j|. One group [14j attempted to explain the slight observable curvature in Fig. 
7(b) by altering the model from s ~ and thus 



This illustrates perfectly why all scientists need to be aware of dimensional issues, for it is impossible 
to take the logarithm of a mass. 'But,' might come the response, 'of course I can. This mouse 
weighs 20g. Look: I type 20 into my calculator, press "log", and there you are'. But a mass has no 
intrinsic number, until we specify a unit. Here this is one gram, and what has actually been done 
is to take the logarithm of the dimensionless ratio of a mass m to a mass unit mo = Ig^ Thus 
the model is actually 



and is no longer invariant under change of unit. If we replace rriQ by rrig, then both a and (3 are 
altered: with fi := log(mo/?TiQ), we have a ^ a + + and /3 1— )• f3-\-2j^. Scale invariance of 
the functional form has been lost, and in particular the linear term, the power /3, is now a function 
of the choice of unit, and can always be set to zero by a correct choice of mo, for (unlike a line) 
a parabola does have a preferred origin. The coefficient 7 of the quadratic term remains invariant, 
and is proportional to the change in gradient of the logarithmic data — but the mere fact of being 
able to fit a non-trivial quadratic is meaningless, for a line segment is approximated arbitrarily well 
by a quadratic, provided its turning point is sufficiently far away (Figure 8) [15]. 

'^Nor does using multiplicative errors solve the problem, for (at a similar mass) we wish to give the same significance 
to a creature with double the expected BMR as to one with half the expected BMR. The logarithmic data do this, 
for I log2| — I log i|, but merely ussing a multiplicative error weighting does not, for 2 — 1 > 1 — |. 

^This explains the precise form of the axis labels in Figs 4, 5, 7 8. In Fig. 4, for example, 'log(i/ft)' means the 
logarithm of the ratio of the length I to one foot. 



log s = a + (3 log m 



(where a and /3 are the constants to be fitted by linear regression) to 



log s = a + (3 log m -|- 7(log m)^ . 
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Figure 8: the quadratic fitted to eutherian data. 

V. The central role of the humble power among the functions used in applied mathematics is 
sometimes masked by the more prominent contributions of the exponential, trigonometric and even 
many of the special functions, the mathematics of which is more substantial. But all science — all 
applied maths — has units and dimensions, and so has dimensionless numbers, and so has scaling 
and thus power laws. 

The initial ideas we have rambled through here lead on to all sorts of profound topics. Power 
law probability distribution functions and the processes which result in them make a lovely story in 
themselves [16]. One of the most mathematically-sophisticated subjects of physics, the theory of 
critical phenomena, has as its output universal classes of such power laws [17]. Most recently, the 
idea that natural processes can enforce criticality (rather than needing to be fine-tuned) has been 
revolutionary [I^. We wish the reader joy in exploring these topics! 
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